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C ^ 3 ' Lattice results on the N to A electromagnetic, axial-vector and pseudoscalar 

form factors are evaluated using dynamical staggered sea quarks and domain 
wall valence quarks for pion masses in the range of 580-350 MeV, as well as, 
(^J ■ dynamical and quenched Wilson fermions for similar pion masses. 
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jy-^ ■ 1. Introduction 

State-of-the-art lattice QCD calculations can yield model independent re- 
sults on N to A transition form factors, thereby providing direct comparison 
with experiment. One such example is the N to A quadrupole form factors 
O |' that have been accurately measured in a series of recent experiments at 

low [1,2] and high momentum transfers [3]. They encode information on 
the deformation of the nucleon and A. We present results on these N to A 
electromagnetic form factors, as well as on the dominant axial- vector N to 
A transition form factors Cf(q 2 ) and C§{q 2 ). Experiments using electro- 
ns) ' production of the A resonance are in the progress [4] to measure the parity 

violating asymmetry in N to A, which, to leading order, is connected to 
Cf(q 2 ). Evaluation of the pseudoscalar 7riVA form factor, G^na^ 2 ), fol- 
lows once the N to A sequential propagators are computed. In addition, we 
evaluate the nucleon axial-vector form factors and the irNN form factor, 
[ GkNn{(1 2 )- Having both the nucleon and the N to A form factors allows us 

to discuss ratios of form factors that are expected to show weaker quark 
mass dependence and be less sensitive to other lattice artifacts. Further- 
more, knowledge of the axial-vector form factors and the ttNN and ttN A 
form factors allows us to check the Goldberger-Treiman relations. 

The light quark regime is studied in two ways: Besides using configura- 
tions with two degenerate flavors of dynamical Wilson fermions we use a 
hybrid combination of domain wall valence quarks, which have chiral sym- 
metry on the lattice, and MILC configurations generated with three flavors 
of staggered sea quarks using the Asqtad improved action [5]. 

2. Lattice Techniques 

Observables in lattice QCD are given by the vacuum expectation value of 
gauge invariant operators in Euclidean time: 

< Q\6\Q >=\\ d[U]d[iP]d[^] 0[U,$, ij;]e- s <>M-s F [u,j,M (1) 



X 

5-H 
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Integrating over the fermionic degrees of freedom we obtain 

<fl\6\n>=^ J d[U] det(Z?[C/])0[ t 7,D- 1 [f/]]e- S9[c/1 (2) 

where D~*[U] substitutes each appearance of —ip n ipj in the operator and 
describes valence quarks whereas det (D[U]) corresponds to sea quarks. The 
path integral over the gauge fields is done numerically by stochastically 
generating a representative ensemble of gauge configurations according to 
the probability pp] = ^ exp { _^ [c/] + ln (det(jD[[/]))} . (3) 

In this work, besides Wilson fermions for the sea and valence quarks, we use 
staggered sea quarks (det(D st aggered[C^])) and domain wall valence fermions 
(Dp W [[/]). The expectation values are obtained by summing over the U- 
ensemble: < Q\6\Q >= limjv->co Z)fcLi 0[U k , D~ 1 [U k ]]. 

The evaluation of form factors involves taking numerically the Fourier 
transform of two- and three-point functions with respect to momentum 
transfer which, on a finite box of spatial length L, takes discrete values in 
units of 2n/L. For large values of momentum transfer the results become 
noisy and therefore we are limited up to Q 2 = —q 2 ~ 2 GeV 2 . To ensure 
that finite volume effects are kept small we take box sizes such that Lm^ ~ 
4.5, where is the pion mass a . In addition, discretization errors due 
to the finite lattice spacing a must be checked. Wilson fermions have 0(a) 
discretization errors and staggered fermions with Asqtad action and domain 
wall fermions (hybrid approach) have 0(a 2 ) errors. Therefore agreement 
between results in these two approaches provides an indication that cut-off 
effects are under control. Finally, we use larger bare u- and d -quark masses 
than physical and extrapolation to the chiral limit must be considered. 

Form factors are extracted from three-point func- 
tions, G AJN (t 2 ,t 1]q ) = < O|E Xl , X2 e lt i Xl f4(x 2 ,i 2 )J(x 1 ,t 1 )jt(0)|r! >, 
shown in the diagram below: 

Interpolating fields for N and A are: 
JP( X ) = e abc [u Ta {x)C l5 d b (x)]u c (x), 




» = ^=e abc {2[u Ta {x)C la d b {x)]u c {x) 
v3 

+ [u Ta {x)C la u b {x)]d c {x)} 
In all cases we apply Gaussian smearing at the source and sink. In the 



a One exception is in the case of dynamical Wilson fermions at the smallest pion mass 
for which Lm w = 3.6 as marked in the Table. 
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case of unqucnchcd Wilson fermions HYP-smearing is applied to the gauge 
links used in the Gaussian smearing of the source and sink. In the case 
of domain wall fermions we use HYP-smeared MILC configurations in all 
computations. We carry out sequential inversions by fixing the quantum 
numbers at the sink and source, which means that the sink time t<i is 
fixed, whereas the insertion time t\ can vary and any operator can be 
inserted at t\. In this work we consider the vector current, = ipj^^tp, 
the axial- vector current, A a ^ = ^7^75 irV' an d the pseudoscalar current, 
P a = ■ipij^^ip, where r a are Pauli matrices acting in flavor space. All 
x\ and x 2 are summed over and we vary t\ in search for a plateau. The 
exponential time dependence and unknown overlaps of the interpolating 
fields with the physical states cancel by dividing the three-point function 
with appropriate combinations of two-point functions [6]. 
The lattice parameters that we use are given in the Table. 



Wilson fermions 



number of confs 


K 


rn w (GeV) 




m N (GeV) 


Quenched 32 3 x 64, = 6.0, a 


- 1 = 2.14(6) GeV 


[a = 0.09 fm) from nucleon mass at chiral limit 


200 




0.1554 


0.563(4) 




1.267(11) 


200 




0.1558 


0.490(4) 




1.190(13) 


200 




0.1562 


0.411(4) 




1.109(13) 






k c =0.1571 


0. 




0.938(9) 




Unquenched [7] 24 3 x 40,3 = 5.6, a 


- 1 = 2.56(10) GeV (a = 0.08 fm) 


185 




0.1575 


0.691(8) 




1.485(18) 


157 




0.1580 


0.509(8) 




1.280(26) 






Unquenched [8] 24 3 x 32,/3 = 


5.6, < 


i" 1 = 2.56(10) GeV 


200 




0.15825 


0.384(8)^- Lrn w 


= 3.6 


1.083(18) 






k c = 0.1585 


0. 




0.938(33) 



Hybrid scheme a 1 = 1.58 GeV (a = 0.125 fm) from MILC collaboration 
number of confs Volume (am„, d ) sea K)° {am q ) DW rn® w (GeV) mjy (GeV) 
150 20 ;i x 64 0.03 0.05 0.0478 0.589(2) 1.392(9) 

198 20 ;i x 64 0.02 0.05 0.0313 0.501(4) 1.255(19) 

100 20 ;i x 64 0.01 0.05 0.0138 0.362(5) 1.138(25) 

150(300 for CMR) 28 ;i x 64 0.01 0.05 0.0138 0.354(2) 1.210(24) 



3. AT to A electromagnetic form factors 

The N to A matrix element of the electromagnetic current can be decom- 
posed into a dominant magnetic dipole, Gmi, and two suppressed electric 
and Coulomb quadrupole form factors, Gei and Gc2- A non-zero Ge2 and 
Gci signal a deformation in the nucleon and/or A. Precise experimental 
data on the quadrupole to dipole ratios, Rem (EMR) = — Q^iiq^ 7 and 

Rsm(CMR) = -^ c"^) , suggest deformation of the nuclcon/A [9]. 
In Fig. 1 we show the EMR and CMR ratios for the smallest pion mass 
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Fig. 1. The EMR (left) and CMR (right) for the lightest pion mass in our three type 
of simulations. 

in the quenched case, for two dynamical flavors of Wilson fermions and 
in the hybrid approach. For the first time in full QCD, we achieve good 
enough accuracy to exclude a zero value for these ratios. Furthermore, at 
low Q 2 , unquenched results become more negative bringing lattice results 
closer to experiment and showing the importance of the pion cloud effects 
at small Q 2 . 



k Hybrid m n =0.50 GeV 
a Hybrid m =0.36 GeV 



O N F =0 m n 

* N F =0 m n 

• N =2 m 



:0.49 GeV 
=0.41 GeV 
=0.51 GeV 



* experiment 



□ N=2 m =0.38 GeV 



0.5 



1.0 

Q 2 



1.5 

(GeV 2 ) 



Fig. 2. 
in the 



Magnetic dipole form factor G* m 
Ash parameterization: G* = 



In Fig. 2 we show lattice results 
for for Wilson fermions and in 
the hybrid approach as well as ex- 
perimental results. At low Q 2 lat- 
tice results are below experiment for 
the pion masses considered here. It 
remains an open question whether 
extrapolation to the physical pion 
mass can bridge the gap between 
lattice and experiment. 



_Q 



4. Nucleon and TV to A axial-vector form factors 

In the case of Wilson fermions, besides N to A we also calculate the nucleon 
axial-vector form factors. The LHP collaboration [10] has evaluated these 
form factors in the hybrid approach with the same parameters as those 
used in our N to A study and therefore, in this case, we use their results to 
compare. The nucleon axial- vector form factors G a and G p are given by 



{N(p')\Al\N(p))=i 



N 



E N (p')E N (p) 



3 

<<■/-') G A (q 2 h^ 5 + 7^G p (q 2 ) T -u{p) 

(4) 

Since the final state is no longer the A a new set of sequential inversions is 
needed. The decomposition of the N to A matrix element of the axial-vector 
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current can be written in terms of four transition form factors [11]: 



m W m A 



A.//-')|A>|^( Pl .)> = V| v ^ 



-A / / fs 
U A+ (P.s) 



+ 6 2 J q\q» + 



m N 



-7 + 



N 



«p(p,s)(5) 



Under the assumptions that ~ and <C C^ 1 the parity violating 
asymmetry is proportional to the ratio C^jC\ [12], where C\ can be 
obtained from the electromagnetic N to A transition. 



N F =0 m n =0.49 GeV 

X N p =0 m n =0 (linear extrap.) 

• N p =2 m[=0.51 GeV 

* N r =2 m n =0 (linear extrap, 

♦ MILC m =0.50 GeV T 

I 



o.o 1 — 

0.0 



i 1 if 




Q 2 (GeV 2 ) 



As shown in Fig. 3, the asymmetry is 
non-zero when Q 2 = [6], increases 
with Q 2 up to about Q 2 ~ 1.5 GeV 2 
and shows small unqucnching effects 
for this range of quark masses. Given 
this weak quark mass dependence, 
the results can be taken as a predic- 
tion for the ratio to be measured by 
the GO collaboration [4]. 



Fig. 3. Parity violating asymmetry. 
4.0.1. Goldberger-Treiman relations 

Partial conservation of axial current (PCAC), d^A® = / T m 2 r 7r a , and the 
axial Ward Identity, d^A a — 2m q P a , relate the pion field 7r a with the 

pseudoscalar density: 7r° = (2m q P a / f^m 2 ) , where the pion decay constant 

fx is determined from the two-point function < 0|^4^|7r b (p) >= ip^8 ab 'f n . 



where 



The renormalized quark mass, m q , is given by m q = "1^1^? wl?"^ ■ 

Aq and P a are the renormalized currents. To obtain the irNN and 7rAA 
form factors we use the decomposition 



2m q (N(p')\P 3 \N(p)) = i 



1 m 2 N 


f^m 2 GnNNiq 2 ) 


E N (p')E N (p) 


ml - q 2 



2m q (A(p , )\P 3 \N(p)} = i^\J 



m A m N hm^G^M^q) 



Ea{v')E N (v) 



ml - q 2 



qu 
2m jy 



u P (p) 



PCAC relates the axial form factors Ga and G p with G^nn and equiv- 
alently and C$ with G^na- These are the well known generalized 
Goldberger-Treiman relations (GTRs). As mentioned above there are ad- 
vantages in considering ratios. In Fig. 4 we show two such ratios, namely 
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GttNa/Gt:NN and 2C^/Ga- Both are independent of Q 2 and the quark 
mass. Fitting to a constant we find 2C^/Ga ~ 1.6 ~ G^na / G^nn , which 
implies the Goldberger-Treiman relations: G^nn^ 2 ) fir = wjvGa^ 2 ) and 
G-nN&iq 2 ) fir — 2mjyC^(q 2 ). Assuming pion-pole dominance we can write 
G P (Q 2 ) = ^f^r G A (Q 2 ) and C£(Q 2 ) = C£(Q 2 ). Therefore 

we have the equality SCq(Q 2 )/G p (Q 2 ) — G^na/G^nn ■ We find that 
8Cq(Q 2 )/G p (Q 2 ) ~ 1.7 [13] a few percentage larger than G^na/G^nn- 

1 
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e 1.5 

| 

1.0 
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' 5 r° N =0 m =0.49 GeV 




-|-f - 
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* N =0 i 
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=0.41 GeV 



A N F =2 GeV 
• N f =2 in=0.51 GeV 

a N =S m =0.38 GeV 



1.0 1.5 

Q 2 (GeV 2 ) 
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1.5 



Q 2 (GeV 2 ) 

Fig. 4. The ratios G^na/G^nn (left) and 2C£/G A (right). 
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1 i 


jr. — 1 


1 1 
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• N T =3 m^O.51 GeV 




n N,=2 rn^O.38 GeV 




* Hybrid m n =0.59 GeV - 




A Hybrid in=0.50 GeV 




H Hybrid m n =0.36 GeV 







1.0 1.5 

Q 2 (GeV 2 ) 



Fig. 5. The nucleoli axial form factors Ga (top, left) and G v (bottom, left) and the N 
to A axial form factors (top, right) and C§ (bottom, right). 

In Fig. 5 we present the nucleon and N to A axial form factors separately 
together with fits of Ga and C§ to a dipole form, go /(-%- + l) 2 . Dynamical 
QCD results in the hybrid approach for the smallest pion mass, where we 
can access low Q 2 values, show large unquenching effects. Having fitted Ga 



February 2, 2008 2:32 



WSPC - Proceedings Trim Size: 9in x 6in JLab'proc2 



7 



and C§ , we can check if pion-pole dominance describes the Q 2 -dependence 
of G p and C$ . The dashed lines correspond to the quenched data and show 
the behavior for G p and C§ extracted from fits to Ga and C§ assuming 
pion-pole dominance, whereas the dotted line shown for C§ is for the hybrid 
approach, in both cases for the lightest pion mass. As can be seen, they 
deviate from the lattice results at low Q 2 . Instead they are best described 
by the solid curves, which are obtained by fitting the pole mass. 



H h 



H h 



N =0 m =0.41 GeV D N =2 m =0.3B GeV 



e I s j 4 

CN 1.0 J" 1 — — 



a Hybrid m =0.36 GeV 




0.2 0.4 0.6 0.8 1.0 1.2 1.4 
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Fig. 6. Left: The ratios Rnn (top) and Rna (bottom) and right: G^nn (top) and 
G^na (bottom) for the smallest pion mass in each type of simulation. 



In Fig. 6 we show the relations 



Rnn 



4mj V / 7r G*j VA f(Q 2 
m 2 G p (Q 2 ) 



Rna 



^Af/7rG'*j VA (Q 2 ) 

2mlC A (Q*) ' 



(6) 



where we have defined G* 7rNN {Q 2 ) = G^NNiQ 2 )/ iX + Q 2 l m V) with a cor- 
responding expression for G* JVA . As can be seen these ratios are con- 
sistent with unity for all Q 2 -values. Finally, in Fig. 6 we show G^nn 
and GkNA for the smallest pion mass. The dash lines are obtained from 
fits of Ga and C$ via the GTRs, G^nnIQ 2 ) = m N GA(Q 2 )/ '/ir and 
GkNa{Q 2 ) = 2mArC^(Q 2 )//7r. As can be seen, there are large devia- 
tions at small Q 2 . Lattice results at this pion mass give a smaller value 
in the limit Q 2 — > than what is extracted from experimental data namely, 

G %NN (0) = 13.21(11) [14]. The solid lines are fits to the form [1 - A^L 
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with a, A fit parameters. Using these fits Gknn(0) and G t atat(0) are ex- 
tracted (for details see Ref. [13]). 

5. Conclusions 

Lattice results on the electromagnetic, axial-vector and pseudoscalar form 
factors for the nucleon and the N to A transition are presented in the 
quenched approximation, for two-flavors of dynamical Wilson fermions and 
using dynamical staggered sea quarks and domain wall valence quarks (hy- 
brid approach) . Results on the quadrupole to dipole ratios EMR and CMR, 
obtained in the hybrid approach reaching down to a pion mass of 350 MeV 
and low Q 2 -values, are non-zero and of similar magnitude as in experiment. 
We also find that ratios of form factors, such as G^na/G^nn ~ 1.6 and 
2C£/G A ~ 1.6, calculated using Wilson fermions, are in agreement with 
phenomenology. Our results for the ratio C^/Cif as a function of Q 2 , can 
be regarded as a lattice prediction for the parity violating asymmetry to 
leading order. The deviations from experiment seen for the magnetic dipole 
N-A transition form factor G* m and the values of G^nn and G^na in the 
limit Q 2 — > need further study. In particular, finite lattice spacing ef- 
fects, as well as, chiral extrapolation to the physical pion mass must be 
investigated. 
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